Andreev spectroscopy and surface density of states for a three-dimensional time-reversal invariant 

topological superconductor 
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A topological superconductor is a fully gapped superconductor that exhibits exotic zero-energy Andreev 
surface states at interfaces with a normal metal. In this paper we investigate the properties of a three-dimensional 
time reversal invariant topological superconductor by means of a two-band model with unconventional pairing 
in both the inter- and intraband channels. Due to the bulk-boundary correspondence the presence of Andreev 
surface states in this system is directly related to the topological structure of the bulk wavefunctions, which 
is characterized by a winding number. Using quasiclassical scattering theory we construct the spectrum of 
the Andreev bound states that appear near the surface and compute the surface density of states for various 
surface orientations. Furthermore, we consider the effects of band splitting, i.e., the breaking of an inversion- 
type symmetry, and demonstrate that in the absence of band splitting there is a direct transition between the 
fully gapped topologically trivial phase and the nontrivial phase, whereas in the presence of band splitting there 
exists a finite region of a gapless nodal superconducting phase between the fully gapped topologically trivial 
and nontrivial phases. 

PACS numbers: 73.43.-f, 73.20.At, 73.20.Fz, 74.25.Fy 



I. INTRODUCTION 

Due to the recent experimental discovery of the quantum 
spin Hall effect^ and the three-dimensional, spin-orbit in- 
duced %2 topological insulator^— there has been a surge of in- 
terest in the study of topological insulating electronic phases. 
Parallel to these developments, many workers have examined 
topological superconductors-12^ which are fully gapped un- 
conventional superconductors that exhibit exotic gapless An- 
dreev surface states. Both topological insulators and topolog- 
ical superconductors can be described within a unified math- 
ematical framework / 'i 15 ' 16 ! 22 which provides a complete and 
exhaustive classification of topological phases of gapped free 
fermion systems in terms of discrete symmetries and spatial 
dimension. A distinctive property of these states is the bulk- 
boundary correspondence, which connects the presence of de- 
localized boundary modes to the topological structure of the 
bulk wavefunctions. 

Notable examples of topological superconductors include 
the spinless chiral (p x +ip y )-w&ve superconductor— and the B 
phase of superfluid 3 Hei U ' 12,17 Here we focus on another type 
of topological superconductor, which has been largely over- 
looked so far, namely the three-dimensional superconductor 
in symmetry class CI in the terminology adopted by Ref. [U 
(following the Altland-Zirnbauer classificatio n 24 ^ 5 ). The dis- 
tinguishing characteristic of the CI topological superconduc- 
tor is that, unlike any of the other three-dimensional topolog- 
ical states, it possesses a form of SU(2) spin- or pseudospin- 
rotation symmetry. In Ref. 14 a tight-binding model on the di- 
amond lattice was proposed that realizes this nontrivial topo- 
logical phase. 

In this paper we recast the model of Ref. 03 into a form 
in which the topology of the Bogoliubov-de Gennes Hamil- 
tonian is completely determined by the phase structure of the 
superconducting gaps near the normal state Fermi surfaces. 
That is, we consider a two-band superconductor with exotic 



inter- and intraband gap functions, whose topological char- 
acteristics do not depend on the full Brillouin zone, but are 
controlled entirely by the properties of both the inter- and in- 
traband gaps in the neighborhood of the Fermi surfaces. The 
reason for considering this case is two-fold: (i) It provides a 
clear interpretation of the topological properties in terms of 
the phase winding of the superconducting gap functions and 
(ii) it allows for the straightforward application of the tools of 
quasiclassical scattering theory, a technique which has proven 
to be extremely useful for the study of the pairing symmetry 
in unconventional superconductors Within this formal- 
ism, both the surface density of states and the spectrum of the 
Andreev bound states can be readily computed. It is known 
that in some unconventional superconductors the presence of 
sub-gap surface bound states leads to zero-energy anomalies 
in the surface density of states^^ We will see that this rule 
also applies to the CI topological superconductor. 

The remainder of the paper is organized as follows. Sec- 
tion describes the model Hamiltonian and its symmetries. 
In Section [III] we introduce a bulk topological invariant and 
compute the phase diagram as a function of band width and 
chemical potential. Section|lV]is concerned with the Andreev 
bound state spectrum and the surface density of states for var- 
ious surface orientations. We conclude with a summary and 
discussion in SectionlVl 



II. MODEL HAMILTONIAN AND SYMMETRIES 



Our starting point is a time-reversal invariant two-band su- 
perconductor on a simple cubic lattice with inter- and intra- 
band pairing, which has the form of a 4 x 4 Bogoliubov- 
de Gennes Hamiltonian. The mean-field Hamiltonian H. = 
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^ fe ^ k H(k)^ft is diagonal in momentum space with 

-$-*>) 

and M/ fe = (afc t ,& fc1 -,aL fci ,&L fc 4.) T , where a k<7 and 6 fc(T de- 
note electron annihilation operators with spin a and momen- 
tum k for band one and two, respectively. The normal state 
Hamiltonian h(k) and the gap matrix S(k) are given by 

«*> " ( 6 o* el) - '« - Jl) ■ <>» 

respectively. Here, the band dispersions are 

Qjk = tj (cos k x + cos k y + cos k z ) — fj,j, ( 1 c) 
with j = 1, 2, the intraband pairing potential is 

A fc = Ad (cosfc x - cosfcj,) + A s , (Id) 
and the interband pairing potential takes the form'' 2 
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FIG. 1: (color online). Phase diagram for the CI topological super- 
conductor, Eq. 0, as a function of bandwidth and chemical potential 
with $o = 0.4, Ad = 0.4, and A s = 0. The gapped phases are char- 
acterized by the even-numbered winding number v, Eq. dlO| >. Blue 
areas are bulk gapless phases, (a) Inversion asymmetric case with 
split bands, t\ =t,t% = 0.9t and fii — fi2- (b) Inversion symmetric 
case with degenerate bands, t\ = ti and \x\ = yb%. 



= $0 (sin k x sin k y + i sin k z ) . (le) 

We observe that Hamiltonian (Q]i is closely related to the 
model of Ref. Specifically, by performing a particle- 
hole transformation the tight-binding Hamiltonian analyzed 
in Ref. 14 can be brought into a form in which the momentum 
dependence of the gap functions along the Fermi surface have 
the same topology as the pairing potentials in Eq. ( flab . We 
note, however, that the present model is defined on a simple 
cubic lattice, whereas in Ref.Qjja diamond lattice was consid- 
ered. The energy spectrum of H(k) is composed of four bands 
with energies E(k) g {-Aifc, -A 2 fe, +A ifc , +A 2 k} and 



Ai fc = ^A 2 + - [9i fc - 9 2fe + B k } 2 , 

A 2fc = JAl + ^[e lk -e 2k -B k } 2 , (2) 

where B k = ^4 |$ fc | 2 + (6 lfc + 9 2fe ) 2 . 

A topological superconductor belonging to the symmetry 
class CI satisfies two independent antiunitary symmetries: 
time-reversal symmetry T = KUt, with T 2 = +1, and 
particle-hole symmetry C = KUc, with C 2 = — 1. Here, 
K denotes the complex conjugation operator. For the above 
example, time-reversal symmetry is expressed as 

U T H*(-k)U^ T = +H{k), (3a) 

where Ut — 1 is the 4 x 4 identity matrix. The particle-hole 
symmetry can be expressed as^ 

U c H*(-k)ul = -H(k), (3b) 

with Uc = i&2 ® Here, and in the following, 0x2,3 de- 
note the three Pauli matrices and er is the 2x2 unit ma- 
trix. By combining Eq. d3ab with Eq. ( l3bb we find that the 



gap matrix S(k) is required to be Hermitian, 5'(k) — S(k). 
Besides particle-hole and time-reversal symmetry, Hamilto- 
nian (Q3 satisfies another symmetry which is given by 

S^H(k)S = -H(k), (3c) 

with S = iUtUc = —<J2 ® 0tj. That is, H(k) anti-commutes 
with the "chiral" symmetry operator S. The significance of 
symmetry d3cb is that it allows us to bring the Hamiltonian 
into block off-diagonal form. Namely, we find that in the basis 
in which S is diagonal, S = WSW^ = diag (—(Jo, +00), the 
Hamiltonian H(k) takes the form 

H{k) = WH(k)Wi = ^ t ° (fe) , (4) 

with the off-diagonal component 

D(k) = h(k) -i6{k), (5) 
and the unitary transformation 

W = 4= h Z*? ) ■ (6) 

The block off-diagonal form © is particularly useful to un- 
cover the topological properties of the ground state wavefunc- 
tions, as we will explain below. 

We end this Section by discussing the case of degenerate 
bands in Eq. ©, i.e., 0(fc) = 0i(fc) = 8 2 (fc). This condi- 
tion leads to an additional symmetry of H(k), 

U p H{-k)Ul = H(k), (7) 

with the unitary matrix U p = 03 <E) cr 2 . Eq. (0 represents 
a type of inversion symmetry, as it is a unitary symmetry 
that relates Bogoliubov-de Gennes Hamiltonians at +k and 



3 



—k via a transformation that interchanges the two bands. 
With symmetry Q the energy spectrum becomes degener- 
ate and takes the simple form E(k) £ {— Afc,+Afc} with 



A, 



III. WINDING NUMBER AND PHASE DIAGRAM 

To determine the topological properties of the model under 
consideration we first introduce an integer topological invari- 
ant, the winding number i/M^ In order to do so it is con- 
venient to adiabatically deform H(k) into a flat-band Hamil- 
tonian. This can be achieved by means of a singular value 
decomposition. First of all we note that for the off-diagonal 
block D(k), Eq. ©, which is in general non-Hermitian, there 
exists a factorization of the form D(k) = U'(k)Ti(k)V(k), 
where S(fc) is a diagonal matrix with positive real numbers on 
the diagonal and U(k) and V(k) are unitary matrices. Direct 
calculation shows that the eigenvalues of S(fe) are identical to 
the positive eigenvalues of the Bogoliubov-de Gennes Hamil- 
tonian H(k). For a fully gapped superconductor, it is possi- 
ble to adiabatically deform the spectrum into flat bands with 
eigenvalues +1 and — 1. This procedure amounts to replac- 
ing by the unit matrix. Hence, the flat-band Hamiltonian 
Q(k) in the off-diagonal basis reads 



Q(k) = 



q(k) 
gt(fc) 



(8) 



with the unitary matrix q(k) = W (k)V(k). In terms of 
the gap functions and band dispersions of model (JTJ, the off- 
diagonal block of the flat-band Hamiltonian is given by 



q(k) 



A+B k - A_6 H 



-2iA_$ fe 
A +J B fc +A_e+ 



D(k)/2 
AifcA 2 fc-Bfc 



(9) 



where A± = An, ± A 2 fc and + = Oifc + 2 fc . As a con- 
sequence of time-reversal invariance q(k) satisfies q T (—k) = 
q(k). The topological invariant characterizing CI topological 
superconductors is defined as the winding number of the off- 
diagonal block g(fe)^iiii 
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d 3 ke^ p Tr [(q-%q) (q^q) (q-%q)] , 



(10) 

where the integral is over the first Brillouin zone. From the 
constraint q T (—k) = q(k) it follows that v is even. 

Next we use the topological invariant ( fTOb to analyze the 
phase diagram of H(k), Eq. (Q~|), as a function of chemical 
potential and band width. Fully gapped phases with different 
topological properties are separated by regions (or lines) of 
nodal superconducting phases (see Fig. |TJ. The condition for 
the existence of a gapless phase can be expressed in terms of a 



vanishing determinant, i.e., det H{k) = 
By use of Eq. ( fTbl and Eq. (f5]l we obtain 



|dct D{k)Y = 0. 



Let us first focus on the inversion asymmetric case with split 
bands, 9ifc 7^ 2 fc. The above two conditions then reduce to 



A fc =0, 



l<I»i 



= -eu,e 



2k- 



(12) 



A 2 k + \<t> k \ + ei fc e 2fc = 0, A fc (Gi* - e 2fc ) = o. (ii) 



Provided sgnOifc = — sgn0 2 fc, Eqs. ( fT2b have solutions 
describing nodal rings that appear in the gapless phase of 
Fig. |TJa) (blue shaded area). These gap-closing lines in mo- 
mentum space are topologically stable and are characterized 
by an integer topological charge, 34 akin to the band touching 
points in graphene. That is, the appearance of these nodal 
lines is generic and stable against (small) perturbations of the 
Hamiltonian, such as, e.g., the inclusion of higher d-wave gap 
harmonics in the intraband pairing potential. 

In the presence of inversion symmetry (O, with Ofc = 
©ifc = ©2fc, the gap closing condition ( fTTT > becomes A| + 
I'I'fcl 2 + ©I = 0. Hence, there are in general four conditions 
that need to be satisfied for the gap to be zero, which exceeds 
the number of free parameters (k x , k y ,k z ). In other words, if 
we consider tuning a single parameter, e.g., the band width t, 
to drive a transition from a topologically nontrivial phase to 
a topologically trivial phase, a gap closing can only occur at 
isolated points in the (k,t) parameter space [see Fig.QJb)]. 
Thus, the presence of inversion symmetry (O leads to a direct 
quantum phase transition between two distinct gapped phases 
(i.e, there is no intervening gapless phase). The phase bound- 
aries in the (i, /i)-plane are given by fx = ±t and /i = ±3i, 
as shown in Fig. 01t>). The fact that direct transitions from 
one gapped phase to another are only possible in the presence 
of inversion symmetry is a feature which is common to three- 
dimensional topological phases. 34 In particular, it also occurs 
in 7L2 topological insulators.— 

To determine the topological nature of the eight gapped 
phases in Fig. [TJ (white and dotted areas), we computed 
the winding number v numerically, by discretizing the inte- 
gral ( [Tol l over the Brillouin zone. It turns out that four phases 
are topologically nontrivial with winding number v = ±2. In 
these nontrivial phases, there appear linearly dispersing, mid- 
gap surface states when the system is placed next to a normal 
metal or insulating state. These exotic Andreev bound states 
are robust against localization from random impurities. We 
will study these surface states in more detail in Section HVl 
In passing, we note that for the inversion symmetric case of 
model ([TJ, there exists an intimate connection between the 
topological properties of the Bogoliubov-de Gennes wave- 
functions, as characterized by the winding number v , and the 
Fermi surface topology in the normal state. In particular, the 
transition from a topologically trivial to a nontrivial phase in 
Fig. |TJb) coincides with a change in Fermi surface topology. 
A similar relationship has been previously reported for fully 
gapped time-reversal invariant spin-triplet superconductors. 18 



A. Fermi surface topological invariant 

An important property of Hamiltonian ([TJ) is that in the pres- 
ence of inversion symmetry its topological characteristics are 
completely determined by the momentum dependence of the 
superconducting gap functions along the normal-state Fermi 
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FIG. 2: (color online). Panel (a): Variation of arg(<£>fc) over the 
Fermi surface. The argument of shows four distinct singularities 
on the Fermi surface. The locations of these singularities are related 
by a four-fold rotational symmetry about the z axis. Panel (b) and 
(c): Variation of sign(Afc) over the Fermi surface. Red indicates 
sign(Afe) = +1, while blue is sign(Afc) = — 1. The ratio between 
s- and d- wave components is A a /Ad — and A s /Ad = 3/7 
in panel (b) and (c), respectively. The normal state band structure 
parameters are ti = t<z = 1 and fj,i = //2 = 1.75. 



surface. To demonstrate this, we give an illustration of the 
topological properties of H(k) in terms of the phase structure 
of the gap functions on the Fermi surface. 

For definitiveness we consider a Fermi surface of spherical 
topology which is centered around the V point, k = 0, i.e. we 
focus on the region t < fj, < 3t in the (i, /i)-plane of Fig.[TJb). 
In the following we will hold the band structure parameters t 
and fi fixed and use the ratio between s- and d-wave compo- 
nents in the intraband pairing ( fTdl ) to tune the system form a 
topologically trivial to a nontrivial phase. Projected onto the 
Fermi surface, we find that $(fe) has four first-order zeroes at 
momenta k±± = (±kp, 0, 0) and &2± = (0, ±fcp, 0), with 
associated singularities in arg($fc) [see Fig. [2ja)]. In other 
words, the real vector field <f>fc over the Fermi surface exhibits 
vortices at ki± and anti-vortices at &2± with winding num- 
ber +1 and —1, respectively. The appearance of vortices in 
<f>(fc) is a necessary but not sufficient condition for the non- 
triviality of the model. What is required in addition is that 
A(fe) reverses sign between vortices of $(fc) with opposite 
winding number. Figs.|2jb) and|2|c) display the variation of 
the sign of A(fe) over the Fermi surface for two different pa- 



rameter choices, which both lead to a nontrivial state. For 
Ad cos(kp) > A s the sign of A& reverses between k\± and 
&2±, an d the system is in the topologically nontrivial phase. 
If A,j cos(fci?) < A s , on the other hand, the sign of A& is the 
same across the Fermi surface, and so we have a topologically 
trivial state. 



IV. ANDREEV BOUND STATES AND SURFACE DENSITY 
OF STATES 

A physical consequence of the nonzero winding number v 
is the appearance of gapless Andreev bound states at the sur- 
face of a CI topological superconductor or at an interface be- 
tween a normal metal and a CI topological superconductor. 
The bulk-boundary correspondence relates the number of An- 
dreev bound states to the topological number v. In this Sec- 
tion we derive the energy spectrum of the bound states and the 
surface density of states using quasiclassical scattering theory. 
For simplicity, we focus on the inversion symmetric case of 
model (Q]i and assume a spherically symmetric Fermi surface. 
But the results we obtain are expected to remain qualitatively 
unchanged upon inclusion of anisotropic Fermi velocities or 
inversion asymmetric perturbations. 

As we have seen in Section ITlTI the topological characteris- 
tics of the inversion symmetric model ([]]) are fully determined 
by the phase structure of the pairing functions on the Fermi 
surface. This implies that we can capture the key topological 
structure of the superconducting state by adopting an effective 
low-energy quasiclassical description as long as the gap func- 
tions A(r) and $(r) are slowly varying over length scales of 
the order of the inverse Fermi momentum kp 1 . Hence, we 
proceed by approximating the momentum dependence of the 
gap functions in the vicinity of the Fermi surface by 

A fc = A[X(kl-k 2 v )/k 2 F + (l-X)}, (13) 
$ fc = t>(k x k y /k F + ik z /k F ) = |$ fc |e^ fc . (14) 

Here we have introduced the parameter A to tune A& from 
d x 2_ y 2-wave symmetry (A = 1, topologically nontrivial) to 
s-wave symmetry (A = 0, topologically trivial). 

We can now apply standard methods to obtain the surface 
bound states. 27 We describe the system in terms of coordinates 
parallel (rii) and normal (a;j_) to the interface. We assume that 
the superconductor occupies the region defined by x± > 0. 
We solve the Andreev equations using the ansatz 

4 

*(k\\,r) = J2 a j*j(k\\,r) (15) 

i=i 

for the wavefunctions of a bound state of energy E, where the 
spinors are written 

*i(fc||,r) = (1, 0, u fcl v k ) T e ik r e- K ^ (16a) 

* 2 (fe||,r) = (1, 0, u £l v- k ) T e fkr e- K ^ (16b) 

* 3 (fc||,r) = (0, 1, v h , ~u k ) T e lk r e-^ x ^{\6c) 

*4(fc||,r) = (0, 1, v k , -u k ) T e ik - r e- K ^(16d) 
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FIG. 3: (color online). Panels (a), (c), and (e): Energy of the upper 
branch of the interface states Ek [see Eq. dl8H as a fraction of the 
gap amplitude for (a) the (100) interface perpendicular to e x , 

(c) the (110) interface perpendicular to -j= (e x + e y ), and (e) the 
(001) interface perpendicular to e z . The color scale is such that blue 
corresponds to E^/Sk — 0, while grey is Ek/S,k = 1- Panels (b), 

(d) , and (f): Comparison of the Andreev state dispersion Ek and the 
bulk gap Hfe for the same interface configurations as in panels (a), 
(c), and (e), respectively. In all panels we take Ad = 1, A s = and 
$o = 0.2. 



the wavevectors k = (kn, k±) and k = (fen, — kj_) with the 
requirement that |fe| = |fe| = kp. 

The energy of the bound states is obtained by the condition 
that the equation \&(feii , J*)Uj_=o = nas a nontrivial solution 
for the coefficients otj. Although in general this yields a rather 
complicated expression for the bound state energies, we can 
simplify matters considerably if we assume that | A h \ = Ag 
and |$fc| = > which holds for certain high-symmetry re- 
flection planes. We hence find the bound state energies 



l r 



~l + A k A~ k + \<S> k \ cos(tp k -<Pk) 



(18) 



This expression is valid in both the topologically trivial and 
nontrivial cases. A zero-energy state is possible whenever the 
following two conditions are both satisfied 

(i) |$fe| = or cos(<£fe - (p k ) = -1, 

(ii) |A fe |=0 or sgnA h sgnA s = -l. (19) 

We shall illustrate different possible combinations of these 
conditions for the appearance of the zero-energy states by ex- 
amining three distinct cases: the (100) surface, the (110) sur- 
face, and the (001) surface. 

The surface bound states in unconventional superconduc- 
tors can be observed by scanning tunneling spectroscopy of 
the surface density of states (SDOS). It is therefore inter- 
esting to consider the SDOS for our three surfaces, in or- 
der to determine the experimentally relevant signatures of 
the topologically nontrivial phase. To obtain the SDOS we 
must first calculate the quasiclassical retarded Green's func- 
tion (r , r '; E). This will not be explicitly constructed here 
as it is rather laborious; for a detailed discussion see Ref.l30l 
and references therein. Because we are dealing with a two- 
band system, the quasiclassical Green's function Q will be a 
4x4 matrix in the Nambu-band space. The SDOS is simply 
the local density of states (DOS) at the surface of the super- 
conductor, where the local DOS at the point r is defined as 

p( E > r ) = -~ E Im {O r ' r; ^ + g h 2(r ' r; *>)} ■ 

fey 

(20) 

Here Q r ' n and Q r ' 22 are the electron-like Jjreen's functions 
for bands one and two, respectively, and E = E + iT con- 
tains the phenomenological broadening parameter T. In all 
our calculations we set T = 0.01 A. 



with 
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— p 
1 



E + i- 



~l-E 2 



(17a) 
(17b) 
(17c) 



and vp = dQp/dp denotes the Fermi velocity. We define 



A. (100) surface 

We consider first the appearance of zero-energy states at 
the (100) surface for A = 1. From the conditions Eq. ( fT9b 
on the intraband pairing we have sgnAfcSgnA^ = 1 for all 
fc, but Afc = along the lines defined by k = (k x , ±k x , k z ). 
The interband potential is always non-zero for these momenta, 
but we do have cos(</3fc — ip-A = —1 (i.e. a sign reversal 
of upon reflection) whenever k z = 0. We hence ob- 
tain two zero-energy states at k\ = (kp/y/2, kp/^/2, 0) and 
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&2 = {kF/V2,—kp/\/2,0). These arguments still hold if we 
decrease A towards 0.5, but the momenta k\ and &2 will move 
towards one another, eventually merging at A = 0.5 where we 
have a gapless state. 

In Fig. [3ja) we plot the A = 1 surface states in units of 
Hfe. The energy of the surface states only deviates signifi- 
cantly from the bulk gap close to the nodal lines of this 
is clearly visible in the k z = cut through the surface states 
shown in Fig.[3jb). 

The SDOS is shown in Fig. Ufa) for fixed A = 1 and various 
values of $. At $ = we recover the SDOS for a 3D d x 2_ y 2 
superconductor. For non-zero $ we observe that at low ener- 
gies the SDOS is linear with a slope that is controlled by the 
velocity of the 2D linearly-dispersing surface states. The min- 
imum energy of the bulk gap is located at E — 0.5$, which 
is visible in the SDOS as the point where the energy depen- 
dence becomes superlinear. The linear slope of the low-energy 
SDOS does not change with increasing $, indicating that the 
velocity of the linear-dispersing surface states is mainly deter- 
mined by A. We note that the finite SDOS at E = is an 
artifact of the broadening parameter. 

As shown in Fig. Ub), the SDOS is qualitatively different 
in the topologically trivial (A < 0.5) and the nontrivial (A > 
0.5) regimes. In the former the low-energy SDOS is vanishing 
within the gap, while in the latter there is the characteristic 
linear energy dependence. Note that the sharp spike at E = A 
in the A = curve is the DOS peak associated with the purely 
s-wave intraband gap. 



B. (110) surface 

We now examine the bound states at the (110) sur- 
face. Again setting A = 1, we see that the condition 
sgn(Afc)sgn(Ajg) = —1 on the intraband potential is satis- 
fied for all fc. A zero-energy state thus only requires that 
vanishes for some k or that cos(ipk — <Pu) = — 1. The former 
holds for ki = (hp, 0, 0) and k^ = (0, kp, 0), while the latter 
is never realized. Compared to the (100) surface the roles of 
Afc and are reversed: the zero-energy states are realized 
because $fc vanishes at momenta where A& has a sign change 
upon reflection. Furthermore, the position of the zero-energy 
states do not change for 1 > A > 0.5, as the sign change of 
Afc at ki and &2 survives while the zeros of remain the 
same. 

The energy of the A = 1 surface states is shown in units of 
the bulk gap in Fig.|3jc). In contrast to the (100) surface, most 
of the interface states have energies differing significantly 
from Hfc. This can be understood as being due to the sign 
reversal upon reflection of A&, which for $fc = would give 
dispersionless zero-energy states for all fci2L2ii In Fig. [2d) we 
show a k z = cut of the interface states and bulk gap. 

In Fig.HJc) we show the change in the SDOS upon varying 
$ at fixed A = 1. Similarly to the (100) surface, at $ = we 
recover the results for a 3D d xy superconductor; note that the 
divergence of the SDOS at E = is due to the zero-energy 
state for all feii. As $ is increased from zero, the SDOS be- 
comes finite at E = and increases linearly with E up to a 
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FIG. 4: (color online) SDOS at the three surfaces as a function of the 
energy E, normalized by the constant SDOS pa of the normal state. 
(100) surface: (a) fixed A = 1 and different values of $ = 0.2nA, 
and (b) fixed $ = 0.5A and different A = 0.2n, offset by 0.25n, 
n = ... 5. (110) surface: (c) fixed A = 1 and different values 
of $ = 0.2nA, and (d) fixed $ = 0.5A and different A = 0.2n, 
offset by 0.75n, n — ... 5. (001) surface: (e) fixed A = 1 and 
different values of $ = 0.2nA, and (f) fixed $ = 0.5A and different 
A = 0.2n, offset by 0.25n, n = . . . 5. 



maximum at the edge of the bulk gap. The low-energy linear 
slope of the SDOS decreases with increasing $, revealing that 
the velocity of the Dirac states increases with <&. This is an- 
ticipated by the result that for this geometry the surface bound 
states appear about the zeroes of 



The change in the SDOS as A is tuned through the topolog- 
ical transition is more subtle than for the (100) surface. As 
can be seen in Fig.Ud), there is relatively little change in the 
low-energy SDOS, with the main feature being that the peak 
at E = 0.5$ becomes sharper as A is decreased. At A < 1 we 
note a drop in the SDOS at E = 0.5A; in the topologically 
nontrivial state, the SDOS is finite on both sides of the drop, 
while in the trivial state it is zero on the higher-energy side. 
That is, the large feature in the SDOS at E < 0.5A in the 
topologically trivial state is due to surface bound states with 
non-zero energy. 
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C. (001) surface 

As there is no sign change of upon reflection from the 
(001) surface, the zero-energy states must be located along the 
nodal lines of A^. The interband pairing potential does 
not vanish for these values of k, but the condition cos(ipk — 
ipz) = —1 is fulfilled when the real part of vanishes, i.e. 
for fc = (k x , 0, k z ) and (0, k y , k z ). For A = 1 this implies a 
single zero-energy bound state at k = (0, 0, kp). As for the 
(100) surface, the origin of this state is due to a zero in A& 
and the sign change of upon reflection. Upon reducing A, 
the single zero-energy state at (0, 0, fci?) splits into two zero- 
energy states at the intersection of the plane (0, k y , k z ) with 
the nodal lines of A&. 

The energies of the surface state deviate most signficantly 
from the bulk along the nodal lines of A&, see Fig.[3je). Un- 
like the other two surfaces, the zero-energy state has quadratic 
dispersion at low energy, as shown by the cut along k x = 
in Fig. Of). For 0.5 < A < 1, however, the two zero-energy 
states have linear dispersion at low energies. 

The SDOS at A = 1 is qualitatively different to the other 
cases because the 2D quadratic surface states contribute a con- 
stant DOS, see Fig. HJe). The height of this constant region 
increases with increasing so that at $ = A the gap is com- 
pletely filled. The edge of the bulk gap at E — 0.5$ is sig- 
nalled by the cusp feature. The constant SDOS within the 
bulk gap is only found at A = 1: as shown in Fig.|4|f), for 
0.5 < A < 1 we find the low-energy linear SDOS character- 
istic of linearly-dispersing zero-energy states. In the topolog- 
ically trivial state the SDOS is vanishing within the gap. 

V. CONCLUSIONS AND OUTLOOK 

In this paper we have discussed the three-dimensional CI 
topological superconductor introduced in Ref. [U We con- 
structed a concrete realization of this topological phase in 
terms of a two-band Bogoliubov-de Gennes Hamiltonian with 
unconventional inter- and intraband pairing potentials. This 
lattice Hamiltonian is just one example of a wider class of 
models that all share the same topological properties. Quite 
generally, one is free to add arbitrary small deformations to 
the Hamiltonian without changing its topological characteris- 
tics, as long as the perturbations do not close the bulk super- 
conducting gap. While we do not know in which specific ma- 
terial the considered tight-binding Hamiltonian could be real- 
ized, it is a convenient canonical model that gives valuable in- 



sight into interesting properties shared by general CI topolog- 
ical superconductors. In the presence of inversion symmetry, 
the topological characteristics of this two-band superconduc- 
tor are fully determined by the momentum dependence of the 
gap functions along the Fermi surface. That is, the topolog- 
ical properties are independent of the electronic band struc- 
ture away from the Fermi surface. We have demonstrated that 
the topological invariant v (winding number) can be related 
to the sign reversal of the intraband gap between vortices in 
the interband gap with opposite winding number. This sim- 
ple criterion could be used in the search for CI topological 
superconductors in real materials. Our results suggest to con- 
sider time-reversal invariant systems with orbital degrees of 
freedom, i.e., multi-band superconductors. 

The CI topological superconductor has exotic Andreev 
bound states at its surface or at an interface with a normal 
metal. These gapless modes are due to the bulk topological 
invariant v, which cannot change as long as the superconduc- 
tor remains fully gapped in the bulk. We have used quasiclas- 
sical scattering theory to study the energy spectrum of these 
Andreev bound states for various surface orientations. An 
important measurement technique to observe Andreev bound 
states in unconventional superconductors is scanning tunnel- 
ing spectroscopy. We therefore computed the surface density 
of states and demonstrated that the presence of Andreev bound 
states leads to pronounced anomalies at low energies in the 
scanning tunneling spectra. These features provide key exper- 
imental signatures of the nontrivial topological character of 
the system. 

Furthermore, it would be interesting to investigate the ef- 
fects of the topological nontriviality of the superconductor on 
other experimental probes, such as tunneling conductance or 
Josephson current. In particular, since our model lends itself 
to the application of quasiclassical scattering techniques, one 
could examine the proximity effects in a junction involving 
a CI topological superconductor and a normal metal or a fer- 
romagnet, for example, or alternatively examine vortex struc- 
tures. We leave these interesting questions for future work. 
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